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Spacecraft  Plasma  M&S 


*  Applications:  Hall  thrusters,  FRC,  LIBS, 

Laser  Plasma  Interaction,  Plasma  Discharges 

*  Complex  physics:  excitation/ionization,  transport, 
radiation,  material,  etc. 


*  Multiple  spatial-temporal  and  density  scales. 


Current  focus: 

Develop  advanced  multiscale 
algorithms  for  plasma  M&S 
in  highly  non-equilibrium 
condition  and  with 
collisional-radiative  kinetics 
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Plasma  M&S 


*  Kinetic  equation:  , 

a,/+v-vi/=-e(/,/)+ecs(/,/) 


Boltzmann  Coll.  Op.:  Q(f,f)  =  If  a{ 


v  —  v* 


R 3  S2 


,co)[f(v')f(v,')-f(v)f(v,)\dcodv' 
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Equilibrium  vdf:  Q(f,  /)  =  0  ->  /(v)  =  ■  exp 

(2;rT  J  27 

—  Fluid  regime:  o' «  1 
—  Kinetic  regime:  £  =  (9(1) 

—  Collisional  plasma:  excitation/ionization,  CE  collisions,  radiation, 
etc. 

*  Methods:  moment  method,  PIC,  DNS. 

*  Challenges: 

-  Multiple  species:  /-►/,  Q(f,/)^YlQ(f»fl) 

—  Dynamical  regime: 
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Collisional-Radiative  (CR)  model 


Non-equilibrium  modeling  of  the  atomic  state 
distribution  function  (ASDF) 

—  Detailed  state-to-state  model  of  atomic  transition, 
i.e.,  excitation,  ionization,  line  radiation,  etc. 

—  Rates  derived  based  on  ab  initio  cross  section. 

Examples:  hypersonic  shocks  in  Ar  &  N2 
Complications: 

—  Accuracy  can  require  many  states 

—  Translational  nonequilibrium 
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Maxwellian  CR 


Hydrogen  model:  En=iH(\-\in2)  in=iH/n 2  g, 

Analytical  rates: 

—  Excitation/deexcitation: 
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—  Line  radiation 

Rate  equation: 
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Level  grouping 


CR  modeling:  level-grouping  ' 

•  Group  effective  rates  of  change 
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•  Internal  structure  of  group  is  assumed  Boltzmann  (Tn) 

—  Piecewise  exponential 

•  This  does  NOT  mean  the  entire  ASDF  is  Boltzmann!! 

•  Group  temperature  must  be  determined  ^additional  conservation 


equation,  e.g..  d£n  .. Afc,^T"  A£,^T*  v*' r ^ 

-dT  =  -N^  \2s2s — z - 2sh'aw)  +  2s2s — 7 - 2sE'&m 

•  DrnroHi  iro9  L">"  *■ 


Procedure? 

—  Solve: 
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with: 


Afn  =  E„,(6-f^,)A,i  =  Ar„ 
d  <A E>„: 


C.(T„) 


iterated 


computed 


AE)„ 


dTn 

— ► 


tabulated 


—  However. . .  -  °(£) 

C.(Tn)  -  o(0 


where 


c  =  e~AE,/Tn 


<H0M0 
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Level  grouping 


CR  modeling:  level-grouping 

•  Other  approaches? 

—  Sub-partitioning:  lowest  level  n0  and  total  J\fn  (no  need  for  Sn) 

T*2  T  sAk)  1  r  -I 

■?«  . . . 


j+7  a/W 

AE).TO  1^-^j  =0(t)/0(£) 

—  Sub-partitioning:  lowest  level  n0  and  upper  distribution  J\fn> 

t*7  r 

XT'*  ^  _ in _  l'r'  n _ „  'll  /rr»\  —  sJr\  t rJ r\ 


o(€)/o(£)  ...fails 


CJ*  ^ 

" "  «75)(A£>.(25)  [N't'9"  n  j 

—  Approximate  Zn  by  expanding  around  mean  energy:  ’XFn  ±jrtf  ng,A/ct. 

-  •(  t  *  •  where  .  A/  37~ 

/"  *'»  J 


sr"/T"  Y1  & 


2n(Tn)  =  €■ 

Kin  L  rn  ~  -n  j 

•  With  n0,  J\fn  partitioning:  1/ln(1+s)  ...fails 

•  With  n0,  Nn.  partitioning:  1/ln(s)  ...succeeds! 

•  Improve  with  successive  iterations...  -,,-v 

K  z;(r„)  =  g'n(ln)e 


A  E’fTn 
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Level  grouping:  Numerical  test 


*  Isothermal  heat  bath: 


10-14  . I 

0  2  4  6  8  10  12  14 

£  (eV) 


Works  very  well  (also  tested  in  cooling  regime). 
Much  better  than  uniform  (standard)  grouping. 
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Level  grouping:  Energy  conservation 


•  Energy  conservation: 

•  Conservation  follows  from  definition  sn  /  r>\; 
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Level  grouping:  Energy  conservation 


Energy  conservation: 

•  Finally... Procedure  shown  to  be  equivalent  to  replacing 
energies  by  “effective”  (condition-dependent)  values  («  EOS) 


Excitation: 


^fLino)  =  (£»o  -  ^no)  *  Ionization:  dfHno)  -  (/„  -  . d(+|(lo) 

^(m'lno)  =  (&m'  ~  '  ^(m'lno) 


**(+!*»*)  —  ***(+[»') 


with 


£^o  —  ft*  ^ 


NOW,  energy  is  conserved 

(down  to  round-off) 


10*  10"r 
time  (see) 
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Multi-fluid  equations 


*  Multi-fluid  model: 

•  5-moment: 


&tPs  +  ^  *  (ps^s)  — 


^t(PsUs)  H~  ^  *  (psUs^s  H~  Ps)  —  Zsens  (E  ~h  us  x  B)  f™ 


dt€s  4-  V  •  (us£s  Ps  •  us)  4-  V  •  qs  —  js  •  E  +  ujg 


ki 


•  Add  Maxwell’s  equations’  r  e 

M  v  x  E  =  -dtB  v  •  E  =  —  (ZiPi  -  nt) 

1 

V  x  B  =  fioi  -  -2^E  V  •  B  =  0 

•  Add  collisions: 

—  Elastic  -  Braginskii  terms 

—  Inelastic  -  Rates  depend  on  both  T  and  relative  velocity 


nnnt  j  I  fnftgcri(g]QiM2)d^idQ2d3vnd3vt 


kj{Te,  |un  -  ue|) 


—  Multi-fluid  CR  model  from  fundamental  principles  being  developed 
(incl.  detailed  balance) 
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Multi-fluid  equations 


*  Electromagnetic  shock:  generalized  Brio-Wu1 

—  FV  with  WENO  reconstruction  and  RK3 


^humlak  &  Loverich,  JCP  2003 
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Two-Fluid  model 


*  Assume  fully  ionized  plasma,  electrostatic  field: 

dtpe  +  V  -  (/?eue)  =  0 

dtpi  +  V  -  ( piUi )  =  0 
5t(/>eue)  +  V  *  (/?eueue  +  j?cI)  = 

dt(piUi)  +  V  -  {piUiUi+pil)  - 
Be  +  V  •  [(£?e  +  Pe)ue]  =  —  V  *  qe  + 

-K3neVeik(Ti 

dtEi  +  V  *  \{Ei  +  Pi)ui] 

*  Gauss’s  law  for  electrostatic  field  v .  e 

*  Transport:  Spitzer  with  flux  limiter 

7 znek2Te 


ik(Ti 


— ( ZiUi-n e) 


q  =  min(f  gFS,q5jj,  qsw) 
I  q,sH  | 

qsw  -  — KeVTe 
Qfs'  =  (f)1^2  kTeVTe 


Iz 


3.22554 (Zi  +  0.24) 
1  +  0.24Z, 
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Electrodynamics 


*  Maxwell’s  equations  for  laser  E  and  B  fields  (different  from  plasma) 


dttE  -  c2V2E  +  c2V(V  •  E)  +  —da  =  0 


1 


€0 


duB  -  crV  B  -  -■ V  x  j  =  0 

£q 


*  Cold  plasma  response  via  Ohm’s  law: 

*  Assuming  monochromatic  wave: 

E(x,t)  =  E(x)e_iwf  +  c.c. 

B(x,  t)  =  B(x)e_lwf  4-  c.c . 
j(x.£)  =  j(x)e_lwt  +  c.c. 

*  Wave  equations  becomes: 

V2E  +  JfcgeE  -  V(V  ■  E)  =  0 
V2B  +  A^sB  +  V  (Ins)  xfvxB)=0 


dtj  +  vj  =  fo^pE 


J 


at 


a  = 

A 


w(l  +  ivjuS) 


Complex  conductivity 
Dielectric  function 
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•  Nonlinear  force  (Hora,  PoF  1985) 

—  General:  fnj  =  j  x  B  +  eoEV  •  E  +  eo  (l  +  V  •  (rf  —  1)EE 

—  Assume  steady-state  and  time-average:  ponderomotive  force 


ip  =  {ini)  =  -V  •  [<T)  -  e0  ( n 2  -  1)  (EE)] 

A 

' -  Maxwell  tensor 


*  Collisional  absorption: 


—  Poynting  theorem:  V  •  S  =  —  j  •  E 

I -  Poynting  vector 


-  Time-average:  ^  =  / j  .  £ 

Assume  1-d  propagation: 


Distribution  A  -  Approved  for  public  release;  Distribution  Unlimited.  PA#  XXXXX 


16 


Laser-Plasma  Interaction 


*  Ion  acceleration  due  to  ponderomotive  force 


x  [/im| 
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Beyond  Maxwellian:  Kn  =  0(1) 


Not-too-far  from  equilibrium  (isotropic) 

—  Discretized  EEDF  yields  rate  equations  for 
discrete  elements  (“bin”) 

—  DB  enforced  at  microscopic  level 

—  High-order,  implicit  and  energy  conserving 

—  More  efficient  compared  to  MCC. 

Far  from  equilibrium 

—  MCC  algorithm  for  inelastic  collisions 

—  Can  resolve  anisotropic  vdf 


r<Tf 


AE 
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drii 

dt 

dfi3 

dt 


j  ... 

...  i  ...  bin  # 

r£i+  Ae 

Ne  /  d£fe(e) 

J£i 

-Ni  fa  fcjjjj)  + 

nu  y 

3 

3 

+n>  y.  - 

—  Drawback:  slow  convergence,  reaction  branching, 
singular  rates,  computational  particle  growth 
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Beyond  Maxwellian:  Particle  Merging 


Scheme  built  on  merging  3+  particles  to  2. 


Mass,  momentum  and  kinetic  energy  are 
exactly  conserved;  Electrostatic  energy  also 
conserved  in  physical  space. 


Split  analogously  defined  by  merging  only 
fractions  of  original  particles. 


Octree  in  velocity  space  inhibits 
thermalization  by  ensuring  only  near 
neighbor  particles  are  merged. 


Higher-moment  conserving  schemes  have 
been  obtained  with  increased  number  of 
merge  result  particles  generated. 
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Probability  Density,  Log10(#  /  (m/s)) 


Particle  Merge:  Importance  of  Octree 
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Probability  Density,  Log10(#/  (m/sf ) 


Particle  Merge:  DC  Breakdown  Case 


Tested  in  3D  ES-PIC  of 
1-KV  DC  Breakdown 

MCC-lonization,  Chain 
Branching  &  Cathode 
Secondary  Emission  causes 
exponential  growth  in  cost 

With  Merging:  Density 
matches  despite  vastly 
different  number  of 
Computational  Particles/Cell 

Negligible  overheard 
demonstrated  in  comparison 
of  wall-clock/iteration  with 
merge  every  iteration 

Enables  direct  control  of 
computational  cost  in  particle 
methods 

Future  Work:  Test  merge  in 
non-Maxwellian  laser  plasma 
test  case 
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Conclusion 


*  Multiscale  algorithms  for  nonequilibrium  flows  with  CR 
kinetics 

—  Level  grouping  schemes  of  electronic  states  of  atoms. 

—  Multi-fluid  equations  developed  to  efficiently  capture 
electron  “hydrodynamics” 

—  Particle  merge/split  for  particle  management,  efficient 
sampling,  inelastic  collisions  ... 

*  Ongoing  works: 

—  Multi-D  simulation  with  level  grouping 
—  Modeling  of  inelastic  collisions  in  multi-fluid 
—  High-order  particle  merging  schemes 
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